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The electromagnetic response of a system with two planes 
per unit cell involves, in addition to the usual intraband con- 
tribution, an added interband term. These transitions affect 
the temperature dependence and the magnitude of the zero 
temperature c-axis penetration depth (A c ). When the inter- 
plane hopping t± is sufficiently small, the interband transi- 
tions dominate the low temperature behaviour of A(T) which 
then does not reflect the linear temperature dependence of 
the intraband term and in comparison becomes quite flat even 
for a d-wave gap. It is in this regime that the pseudogap was 
found in our previous normal state calculations of the c axis 
conductivity, and the effects are connected. 

74.25.Nf,74.25.Jb,74.72.-h 



I. INTRODUCTION 

Studies of the in plane low temperature penetra- 
tion depth of single crystals of twinned and untwinned 
YBa2Cu30a; (YBCCv) at optimum doping (a; r^-,6.93) 
have revealed a linear temperature dependence.tTEl This 
has been confirmed in more *ecent measurements on sin- 
gle crystals byJVlao et al.El and on thin films by de 
Vaulchier et al.Q although reports of other power laws 
persist B Extensions to underdoped samples (YBC06.6) 
as well as to the overdoped case YBCO6.99 confirm that 
this linear law does not depend on oxygen doping al- 
though the value of its slope is affected. In untwinned 
single crystals, it has also been found that the zero 
temperature value of the in-plane penetration depth is 
highly anisotropic: A& = 1600 Awhile A a = 1030Awhere 
b denotes the chain direction, and the sample doping 
is near optimal (x = 6.95 with T c = 93.2 K). This 
large anisotropy in X a b(T = 0) is, however, not reflected 
strongly in the temperature dependence of the normal- 
ized penetration depth Xi(0) / XiXX-which is found to be 
nearly the same for i — a and 6jjaE9 with that along the 
chain showing a slightly steeper slope. 

A linear temperature behaviour is, of course, what is 
expected for a superconductor with a gap which has d- 
wave symmetry with zeros crossing the Fermi surface. 
Also, the addition of Zn impuritiesjj which are believed to 
act as unitary scatters, is found to alter the temperature 
dependence of X a b from linear to T 2 , which is expected 
for d-wave. There are now many other experiments that 
confirm the <i-wave symmetry pfithf gap- Examples are 
angular resolved photoemissiorJl2riI3 (ARPES) as well as 



SQUID&EI experiments. These latter experiments have 
been devised specifically to measure the phase of the gap 
rather than simply probe for zeros crossing the Fermi 
surface as does ARPES and penetration depth. 

Detailed theoretical studies of the a-b plane aaisotropy 
have been carried out by the present authoralaEJ in a 
two plane model, coupled through a transverse tunnelling 
matrix element (t±). The general conclusion of these 
studies is that there is a significant amount of conden- 
sate on the chains and that the gap on the CuO chains 
must be of the same order of magnitude as on the Cu02 
planes. The condensate on the chains is, of course, the 
ultimate source of the orthorhombicity which is reflected 
in the observed large in plane anisotropy for A a and A& 
at T=0. r^he recent theoretical work of QIDonovan and 
Carbottecil and of Xiang and WheatJej£3 further con- 
firm these general conclusions. CITScj'EJ experiments in 
which the gap is measured separately on the chains and 
on the plane give added evidence that both gap values 
are of roughly of the same size. An anisotropy-pf-simi- 
lar magnitude is alsp_seen in the DC resistivityE3l23 and 
optical conductivity^, r, 

Recently, Hardy et al.,u in addition to summarizing 
their result for a-b plane anisotropy in the underdoped, 
optimally doped and -Qverdoped cases, have also pre- 
sented c-axis dataJaM'El For YBC0 6 . 6 (T c — 59 K), they 
estimate A c /A a - 31. For YBC0 6 . 95 (T c = 93.2 K) 
A c /A a ~ 7 which is about the same as for the overdoped 
case with x=6.99 (T =89K). They give detailed results 
for the temperature variation of the ratio A c (0)/A c (T) 
and find, in all cases, that at low T, the curve for the 
c-axis falls well above the linear dependence found for A a 
and Xb and may vary as a higher than linear power. A 
striking fact is that the underdoped and overdoped case 
are found to have nearly the same dependence on reduced 
temperature T/T c while the optimally doped case falls 
slightly below. A possible explanation for the rather flat 
behaviour observed in the underdoped case is that it is 
representative of incoherent coupling between the chains 
and planes. This view is consistent with decreasing cou- 
pling between planar elements as oxygen is removed from 
the sample. 

For fully oxygenated YBCO^, the experimental results 
on the c-axis DC resistivity now appear to converge. It 
is found that p(T) is metallic in the c-dkeptjon and has 
a nearly linear temperature dependence. E^lO'El Further, 
the anisotropy ratio p c /pab is of order 30-70 which is con- 
sistent with the value of A c / A a quoted above if we assume 
that an effective mass model is qualitatively valid. On 
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the other hand, when the sample is deoxygenated, the 
c-axis behaviour can change radically with the DC resis- 
tivity p c {T) showing a semiconducting-like behaviour at 
low temperature (dp c /dT < 0). n n 

Also, the ratio p c /p a b is much bigger&til in the under- 
doped samples thas-for optimal doping. For example, at 
T = 100 K, p c /pab^ is of order 10 3 . This large value is 
in accordance with the value found for the corresponding 
ratio of (A c /Ab) 2 . It is to be noted that as the temper- 
ature is increased, metallic behaviour is recovered even 
for YBC0 6 .6- The very large value of p c observed in the 
underdoped case and its semiconducting behaviour has 
been taken as evidence that the c-axis transport can not 
be described by conventional three dimensional coherent 
Bloch transport and that it is dominated by a different 
(incoherent) mechanism due to the highly two dimen- 
sional nature of the electronic states. Evidence for this is 
provided by estimates of the c-axis mean free path which 
find that it is of the order or less-tkan the c-axis spacing. 

Various possible mechanismsEjES for c-axis transport 
have been reviewed by Zha, Cooper and Pines £3 In 
particular, Rojo and LevinEa have considered explicitly 
static and dynamic off diagonal disorders as well as the 
effect of the interplane matrix element i j_ and also include 
several limiting cases. In particular, ordinary transport 
and the pure incoherent case with t± — and static 
impurity scattering are treated. This latter case has 
been discussed in detail in the work of Graf, Rainer and 
Sauls.ta It leads to a c-axis resistivity that is inversely 
related to its in-plane value. The source of incoherent 
scattering in such models is, however, not clear. An- 
other interesting idea is the .dynamic thermal dephasing 
process suggested by LeggettEa in which in-plane thermal 
fluctuations can be larger than t± and so destroy coherent 
c-axis transport. This mechanism, however, cannot ex- 
plain the observation that metallic transport is restored 
at high temperature for the underdoped samples. A fur- 
ther model based on resonant tunneling through a barrier 
layer which due to Abrikosov,c3 models the c-axis d.c. 
resistivity with some succes. Yet another approach as- 
sumes that the copper oxide planes cannot be described 
by Fermi liquid theory-and spin and charge degrees of 
freedom are separated. EJ No consensus has yet emerged 
about c-axis coupling. 

In this paper, we will be mainly interested in the c- 
axis penetration depth. In this case, Radtke, Kostur 
and Levinc3 have already given results for the T = 
value and temperature dependence of A C (T) in a model 
which contains three distinct contributions dependent on 
t± , a pure incoherent impurity assisted contribution as 
in the work of Graf et alE3 and a phonon assisted in- 
terlayer inelastic part. They find that when the direct 
contribution dominates, A c and A a & will have the same 
linear temperature dependence at low T while for im- 
purity assisted hopping Aj?(T) will show a much flatter, 
more s-wave low temperature behaviour as is observed 
in experiments. There are two problems with this sce- 
nario as applied to underdoped YBCO. First, while t± 



may be small, the coherent part may not be negligible. 
Estimates of t± provided by Zha et al.EI give t± ~ 3.0 
meV for YBC06.68- This value is reduced by an order 
of magnitude from t± ~30-40 meV estimated in opti- 
mally doped YBCO but is still significant and is cer- 
tainly much larger than the value of 0.1 meV estimated 
for Bi2Sr2CaCu20s. Secondly, the impurity assisted con- 
tribution vanishes by symmetry for a cZ-wave supercon- 
ductor if the impurity potential is isotropic in momen- 
tum space, i.e. a J-function in direct space. This may 
not be a serious obstacle, particularly in YBCO which is 
orthorhombic and so cannot have a purely d-wave order 
parameter. A more serious problem is the microscopic 
origin of this impurity scattering. This is particularly 
problematic since many experiments would lead us to 
believe that these crystals are in the clean limit as far as 
ab plane properties are concerned. 

In this paper, we study a different mechanism for the c- 
axis penetration depth. We consider the possibility that 
it is a manifestation of the interband transitions which 
become more important in the c-axis conductivity as t± 
is reduced in magnitude. Iri|-»revious work on the nor- 
mal state A.C. conductivity]! 3 ] we found that the usual 
Drude-like intraband contribution is of order t\ and the 
interband contribution goes like tj_ for the in-plane con- 
ductivity. On the other hand, for the c-axis conductivity, 
the intraband contribution goes like tj_ while the inter- 
band remains of order t\. It is clear then that, for suf- 
ficiently small values of t± , the interband contribution 
to the penetration depth will dominate, and A C (T) may 
be significantly different from A a &. It is precisely for this 
case that we found a pseudo gap in the normal state c 
axis conductivity. Thus the apperence of a pseudogap 
and a fiat low temperature dependence in the c-axis pen- 
etration depth are connected in our theory. 

In section ||, we present the necessary formalism. Nu- 
merical results are given and discussed in section III. A 
brief conclusion is included in section pVj . 



II. DERIVATIONS 

To remain simple, we consider two isolated planes cou- 
pled with a transverse matrix element t± . The uncoupled 
planes have dispersion 

£i = —2a\[cos(k x ) + cos(k y ) — 2B cos(k x ) cos(fc y )] — p\, 

(1) 

representing the tetragonal Cu02 plane and 

6 = -2(72 cos(ky) - p 2 (2) 

representing the orthorhombic CuO chains. In (|l|) and 
(Q), 8<7i and 4<72 are the band widths, B is a next near- 
est neighbour hopping amplitude in units o\, and pi and 
P2 are chemical potentials related to the fillings on the 



2 



chains and planes. This model is already quite compli- 
cated to deal with and final results are obtained only 
numerically. Nevertheless, it needs to be recognized that 
it is still oversimplified and may not apply directly to 
YBCO. For instance, in YBCO, there are three planes 
per unit cell: a bilayer of CuC>2 planes and a layer of 
CuO chains in between them. These 3 planes form the 
basic unit cell which is then coupled to the other unit 
cells through additional, presumably insulating, layers. 
Under such circumstances, estimates of the transverse 
coupling t± provided from DC resistivity data may well 
reflect the value of the interccll coupling (which could be 
incoherent) rather than the plane-chain coupling within 
the cell (which could be coherent). Treating such com- 
plications is, however, beyond the scope of this paper. 
Here, we will assume only 2 planes but take one to be 
chains so as to introduce the observed orthorhombicity. 
The transverse coupling ij_ therefore represents both the 
chain-plane coupling and the intercell coupling. 
The mean field Hamiltonian for our model is 



(3) 



a four by four matrix related to the bare vertex function 
7 M (k + q, k) through the unitary transformation U(k) 
which diagonalizes the Hamiltonian. That is 



7a* 



(k, k + q) = U\k) lu (k, k + q)U(k + q). 



(6) 



The quasiparticle energies Ei(k) are the eigenvalues of 
h(k): E 1 = E+, E 2 = -E+, E 3 = E A = 
where 



E± = V4 + A k 
and e± are the normal state band energies 

6 + 6 



± 



6-6 



t 2 . 



(7) 



(8) 



For a static external electromagnetic vector potential 
A(r), the current is given by 



J M (r) = A, (r) 



(9) 



where C k = [cJ kT , ci_ kJ ,, c\ kV c 2 _kj.] and cj kCT creates an 
electron of spin a and momentum k in the plane (i = 1) 
or chain (i = 2) layers. The Hamiltonian matrix is 



h(k) = 



6 Ak i(k) 

A k -6 -t(k) 

t(k) £2 A k 

-t(k) A k -6 



(4) 



In Eq. (f|) t(k) = -2t±cos(k z d/2), which follows from 
a tight-binding model of c-axis coupling. The unit cell 
size is d along the c-axis. The superconducting gap Ak 
is chosen to be the same in both the plane and chain lay- 
ers. We will not determine Ak through solution of a gap 
equation but will instead be guided by experiments which 
have indicated that the gap on the planes and chains are 
very similar. For our simple model we have a single d- 
wave gap characteristic of the entire system. We make 
no attempt to relate this gap to the basic microscopic 
interactions but simply postulate its existence. 

In a previous paper, we derived a general expression 
for the penetration depth of a system with two atoms 
per unit cell [(-where the two atoms form the chain and 
plane layers). E2I The expression for the penetration depth 
required the two particle Green's function: 



where (//, v) indicates spatial directions and the kernel 
K^y is given by 



-<W t^( k ' k + q )]y ft/* ( k + q> k )b* 



/[^(^-/[^(k + q)] 
Ei(k) -Tluj - E^k + q)' 



(5) 



where f(x) is the Fermi distribution function, e is the 
electronic charge and the sum over i and j range from 1 
to 4. In (||), the electromagnetic vertex 7 M (k + q, k) is 



K„ 



1 



(G^(0,0,0)| A=o -G^(0,0,0)} 



(10) 



In equation (|10|), c is the velocity of light, ft is the volume 
and the first term in the curly bracket is the same as the 
second except that it is to be evaluated in the normal 
state with gap set equal to zero. Note that only the 
q = k = and uj = limit has been taken in (|o|) . This is 
sufficient for calculations of the penetration depth which 
follows from 
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(11) 



Our fundamental formula therefore, has the form 

,,2 



G w (0,0,0) = --^[ 7M (k,k)] 

dm) 

J BE; [i ~ 



E.; — Ea 



(12) 



The matrix ?7(k) which diagonalizes the Hamiltonian 
matrix can be found using a two step process. We first 
introduce a unitary matrix Un which would diagonalize 
the normal state alone. It has the form 



U N = 








v„ 


■ 





Un 





V n 







-U n 








Vn 





-U n . 



(13) 



with 



3 



and 



e+-6 



(14) 



(15) 



Then 



ti = ulh U N 



N 

e+ A k 






e_ 
A k 






A k 

— e_ 



Matrix (|T^) can finally be diagonalized by application 
of a second unitary transformation Us to get 



Ei = 



Ul ti U s 



(17) 



where 



U S = 



u+ 


v + 





o ■ 


v+ 


— u + 














U- 










V- 


— 7i_ 



(18) 



with 



and 



! E±+e± 
2E ± 



< E± - e± 
2£± 



(19) 



(20) 



where we stress again, a single gap has been introduced 
for the condensate on both the CuC>2 planes and CuO 
chains equally. This assumption was justified in our pre- 
vious workE3 and has simplified the mathematics enor- 
mously here. 

There remains only to compute the transformed elec- 
tromagnetic vertex 7 = UgUj^jU jvfej which appears in 
equation (E2j). The bare vertices areE2l 



v lx 

di, 







v ly 

v ly 

V 2 y 

v 2y 



(21) 



(22) 



with 



v± 

v ± 

v± 

v± 



1 



and 



4,1 h dk„ 

1 9t 

^ h dk z 

(16) It is straightforward to show that 



(23) 



(24) 



(25) 



w M+ ai^ 







where 



and 



- 1 de± 



-Vlx, 



-{Vly+V 2 y), 



6-6 



(26) 

(27) 

(28a) 
(28b) 
(28c) 



e + - e 

The coherence factors a and /3 have the simple form 
9 1 E+E- + e+e_ + A 2 



0- 



E+E- 



2 _ 1 E+E- - e+e_ - A 2 



(29a) 
(29b) 



2 

Now that 7 has been specified, may be evaluated. 

III. RESULTS AND DISCUSSION 

We start by rewriting equation (jl2j) (which gives the 
penetration depth X(T) through ( |l0| ) and (pi)) in the 
form 



g/Cgi) 



Eft 



2 f{Ei)-f(Ej) 
^ E, - E< 



(30) 



4 



The first term in (30) is the very familiar intraband con- 
tribution and can be rewritten to read 




dk u 



g/Cgj 
dE+ 



dk u 



g/Cg-) 



(31) 



which is just the usual expression for the penetration 
depth in a 2 band model with energies given by equation 
(||). For t± = equation (|3l]) will make no contribution 
to A 2 and we recover the result for 2 decoupled bands. 
It is to be noted that the second term in (j30|), which is 
new and which gives the interband transitions, will also 
give no contribution to A 2 when t± = as it must. This 
can be seen by noting that each of the vertex elements 
with i ^ j, is proportional to t±. The interband 
contribution to coming from the second term in (^) can 
be written as 



CI 



,J(E + )-f(E_) 



E 4 



E_ 



-0- 



2 1 -/(£+) -/(£_) 



E, 



(32) 



where a and (3 are given by Eqs. (29a 



and ( pb[ ). 

These terms describe interband transitions. The term 
with [f(E + ) - f(E-)]/{E + - £L) vanishes at T = 
and increases with increasing T while the one with 
[1 - f(E+) - /(£_)]/(£+ + E-) is 1 at T = and 
decreases with increasing T. This competition between 
these two terms leads to a much flatter temperature de- 
pendence for A _2 (T) than does the more common intra- 
band term (|3l|). Choosing the z direction, we note tha t 
the intraband part ( |3l| ) goes like t 4 because v z in (28c) 
goes like t\ while the interband contribution (32) goes 
only like t \ . For sufficiently small values of t± , it is clear 
that the interband contribution will dominate the c-axis 
penetration depth and that our intuition based only on 
the intraband band term will no longer hold. Our numer- 
ical work, to be described below, will tell us the value of 
t± where this new regime of behaviour is reached. Be- 
fore presenting results it is interesting to note that the 
in-plane situation is quite different with regards to the 
dominance, or lack there of, of the interband processes. 
This can be seen from equation (J3TJ) and the realization 
that, for the x or y directions, the intraband contribu- 
tion is, to leading order, independent of t±, while from 
(p2|), we see that the interband part is proportional to 
Thus, for small t±, the interband contribution to the 
in-plane penetration depth will not be important com- 
pared with the intraband contribution and the familiar 
behaviour will result. On the other hand, when t± gets 
large, the interband contribution to the in-plane pene- 
tration depth can become of some importance. 

In Fig. H, we show our model band structure. We have 
taken {cr u a 2 , Mi, M2, t±} = {70, 100, -65, -175, 20} meV 
with B = 0.45. The dashed line gives the plane and chain 
Fermi surface contour in the 2-dimensional Cu02 plane 
Brillouin zone as a function of momentum components 




FIG. 1. Fermi surface for the chain-plane model. The pa- 
ram- 
eters are {cri, (72, fi\, /Lt2, t±} = {70,100,-65,-175,20} meV 
and B — 0.45. The dashed lines apply to the unperturbed 
case t± = and represent uncoupled plane and chain layers. 
They also describe the case k z = ir/d in equation (^|), i.e. 
the Brillouin zone boundary. The solid curves are for k z = 
in (^|), i.e. the central plane in the B.Z. while the area be- 
tween solid and dashed curves represent the dispersion in the 
z-direction. 



k x and k y . When the transverse matrix element t± is 
taken into account, the Fermi surface acquires dispersion 
in the z-direction. When k z is zero, t(k z ) = 2t and we 
get the hybridized Fermi surface shown as the solid con- 
tours in Fig. |l]. On the other hand, for k z at the zone 
edge (k z — ir/d), t(k z ) = and we recover the dashed 
contours. For other values of k z , we obtain Fermi surface 
contours intermediate between solid and dashed curve 
and the area between these 2 contours gives the z disper- 
sion of our model. This simplified Fermi surface will be 
kept fixed in our calculations and we present results for 
the penetration depth as a function of temperature for 
all three directions x, y, z (or, equivalently a, b, c). 

What is given in Figs, || and || are the temperature 
dependences of [A AI (0)/A p (T)] 2 for t ± = 20 meV, 10 meV 
and 5 meV, respectively. In all cases, we have used a 
d-wave gap of the form 



A k = A [cos(fc x ) - cos(fey)] 



(33) 



over the 2-dimensional Brillouin zone with Aq(T) taken 
to have a BCS temperature variation but scaled so that 
2A (T = 0) /ksTc = 7. This value is chosen to agree with 
experiments on YBCO which find a rather large value for 
the gap to critical temperature ratio. This is also neces- 
sary in order to get an in-plane temperature dependence 
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FIG. 2. The inverse square penetration depth normalized 
to its value at T = in the a (dashed curve), b (dot-dashed 
curve) and c (solid curve) directions as a function of reduced 
temperature T/T c . The transverse matrix element coupling 
the chain and plane layers is t± — 20 meV. 



for the penetration depth which is in reasonable accord 
with experiment. Before describing our results and com- 
menting on them, we stress that in YBCO x as a function 
of oxygen doping, the Fermi surface of the chains will 
certainly be changed in important ways because oxygen 
depletion occurs on the chains as we go from the over- 
doped case x = 7 to the underdoped case with x = 6.6. 
At the same time, the number of holes on the planes as 
well as on the chains will change so the value of the chem- 
ical potential and displacement between bands, \i\ and 
fi2 will change as will the filling factor in each band. Such 
details are not yet known and in the absence of such infor- 
mation, we have chosen to keep the bands fixed as well as 
the filling which was not chosen to specifically model the 
overdoped or underdoped case. Here, we are primarily 
interested in understanding the effect of interband tran- 
sitions on the penetration depth and their variation with 
chain-plane coupling t± which may change significantly 
with doping although again this is not well known. 

In Figs. Q to ||, the dashed curve is for the a-direction, 
the dot-dashed curve is for the 6-direction and the solid 
for the c-direction. In the first case, t± = 20 meV and all 
three directions show about the same temperature depen- 
dence with A^~ 2 having the steepest slope at low T, A~ 2 
having a slightly flatter slope and X~ 2 having the flattest 
slope. This order is as observed in experiment. As the 
transverse hopping matrix element t± is reduced in size, 
we see that the in plane behaviour is not strongly affected 
but the c-axis is. It becomes flatter at low temperature 
and begins to look much more like a usual s-wave case. 



0.2 



0.4 0.6 
T/T 



0.8 



FIG. 3. As in Fig. |2j, penetration depth in the a (dashed 
curve), b (dot-dashed curve) and c (solid curve) directions 
but with tj_ — 10 meV. 



This different behaviour results when the interband terms 
([32|) start becoming dominant over the more conventional 
intraband term (pl|). We might expect this to occur in 
the underdoped situation in YBCO where t±_ is expected 
to be of order 5 meV or so. For optimally doped YBCO, 
interband terms contribute less to the c-axis behaviour 
and the system becomes more three dimensional. 

In Fig. ||, we show the absolute value of the inverse 
square of the zero temperature penetration depth A(0) -2 
for the three directions. Note the difference in scale be- 
tween \~ 2 and A~ h 2 . In the plot of A c (0), we can see very 
clearly the t\ dependence discussed earlier. The magni- 
tude of A~ 2 is smaller than A^ 2 , roughly by a factor of 
two because of the additional screening current carried 
by the chains. The difference is reduced as t± increases. 
This is attributed to the increase significance of the inter- 
band contribution to the 6-direction penetration depth as 
t± is increased. In the x direction, the intraband term in 
(|3l| ) goes like 



l + O 



t- 



for most of the Brillouin zone and the second term is 
small as is the interband term ( |32"|) which also varies as 
0(t 2 /(£i —£ y 2) 2 )- In the region where the plane and chain 
Fermi surface crosses this argument does not apply, but 
v\ x is small there and so we can conclude that the terms 
in t± are not important. This does not hold, however, for 
the y-direction because in that case v\ y and V2 y are not 
necessarily small in regions of the Brillouin zone where 
the corrections in t±_ arc significant and we expect some 
contribution to A~ 2 (T). 
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0.2 0.4 0.6 0.8 1 



T/T c 

FIG. 4. As in Fig. penetration depth in the a (dashed 
curve), b (dot-dashed curve) and c (solid curve) directions 
but with t± = 5 meV. 



IV. CONCLUSION 

We have considered a model of two conducting planes 
per unit cell, one with tetragonal symmetry representing 
a CuC>2 plane, the other with orthorhombic symmetry, 
representing a CuO chain. A transverse hopping matrix 
element t± is included which couples the two bands co- 
herently and introduces the third dimension. For such a 
model, new transitions enter the electromagnetic proper- 
ties that would not be present in the case of one atom 
per unit cell. In addition to the usual intraband terms, 
there are terms which now correspond to interband tran- 
sitions. We have investigated the effect of these terms 
on the value of the zero temperature penetration depth 
and on its temperature dependence. It is found that for 
the a-direction, the effects of t± are small while for the 
6-direction, along the chains, there can be significant cor- 
rections to the intraband terms, although to leading order 
Xb is still independent of t±. In this case, the corrections 
for the interband terms go like t\. 

For the c-direction (perpendicular to the planes), the 
situation is quite different. The conventional intraband 
terms go like t\ while the interband ones go like t\ so, 
for sufficiently small values of t±, the interband contri- 
bution eventually must take over. We find that these 
have a new temperature dependence quite distinct from 
the linear law expected for a c£-wave superconductor at 
low temperature. They give a much flatter curve at low 
temperature for the c-axis superfluid density. Numeri- 
cal calculations indicate that the interband term will be 
important when t± is of order 5-10 meV. This value is 
reasonable for underdoped YBCO. In our previous work 




t ± (meV) 

FIG. 5. Magnitude of the inverse penetration depth at 
T — as a function of the plane-chain coupling. The 
three curves are A a (0) (dashed), A;,(0) (dot-dashed) and A c (0) 
(solid). 



on the normal state, we found that a pseudogap develops 
in the c axis optical conductivity for such small values 
of t±. This offers an alternate interpretation of the c- 
axis data in the underdoped case, which does not rely on 
any incoherent elastic scattering or c-axis out-of-plane 
phonon assisted process. Our calculations naturally give 
a linear dependence in a- and b- directions and a much 
flatter dependence for the c-axis as is observed. On the 
other hand, as the coupling between the two planes is 
increased, the system becomes much more 3 dimensional 
and the temperature dependence of all three penetration 
depths become similar and are linear at low T for a d- 
wave gap. These findings are all in general qualitative 
agreement with the data in YBCO. 

As discussed in the text, however, our results can prob- 
ably not be applied to the oxides without modifications 
because our electronic structure model is too simple. For 
example, the unit cell in YBCO contains a Cu02 bilayer 
as well as CuO chains. The coupling between planes 
within a cell is likely to be quite different from that be- 
tween two cells while in our simplified model a single t ± 
enters. Still, the model allows us to introduce and under- 
stand better the role that interband transitions can play 
in the c-axis penetration depth. It exhibits many of the 
qualitative properties observed in YBCO. 
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